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Cosmic gamma-ray bursts with primary-photon energies > 10 GeV are sought in the data from the 
Andyrchy array obtained in the mode of detection of a single cosmic-ray component during a net observation 
period of 2005.4 d. The distribution of fluctuations of the detector counting rate agrees with the expected 
cosmic-ray background, the only exception being an event with a deviation of 7.9a. Constraints on the number 
density of evaporating primordial black holes in a local region of the Galaxy are obtained for the chromospheric 
evaporation models. 
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Primordial black holes (PBHs) can be formed in 
5j [the early Universe due to the gravitational collapse of 
primordial cosmological density fluctuations. Theoreti- 
cal predictions of the probability of PBH formation are 
strongly dependent on the accepted theory of gravita- 
tion and the model of gravitational collapse. The pro- 
'cess of black-hole evaporation, which is used to seek 
them experimentally, is also far from being completely 
understood. Thus, PBH detection will provide impor- 
tant information on the early Universe and can be a 
unique test of the general theory of relativity, cosmol- 
ogy, and quantum gravity [T]. Knowledge of the PBH 
spatial distribution is important for their direct search. 
J^j 'The number density of PBHs in our Galaxy can be many 
.orders of magnitude higher than the average number 
03 density of PBHs in the Universe [2], so that the con- 
straints based on the results of direct searches can be 
much more stringent than those implied by measuring 
the diffuse extra-Galactic 7-ray background. 

The bursts of high-energy 7 rays are generated at the 
final stages of PBH evaporation. Since the calculated 
temporal and energy characteristics of such bursts de- 
pend on the accepted theoretical evaporation model [3] , 
the experimental technique used to detect PBH emis- 
sion and, correspondingly, the resulting constraints on 
the number density of PBHs in the local Universe are 
strongly model dependent. In the model without a chro- 
mosphere [4], the evaporating particles do not inter- 
act, all quarks propagate freely, and evaporating par- 
ticles fragment independently of each other. The pho- 
ton spectrum is formed by the quark fragmentation and 
decay of unstable hadrons; hence, this spectrum is not 



thermal. According to the DK02 5J and H97 [6J chro- 
mospheric models, the interacting evaporated particles 
form a quasi-chromosphere, which results in the strong 
energy fragmentation and, accordingly, in the steep pho- 
ton spectrum at high energies. The spectra of 7 rays 
emitted by PBHs depend on the time until the end of the 
black-hole evaporation. The evaporation models were 
analyzed in more detail in [3j. The time interval until 
the end of black-hole evaporation during which 99% of 
the 7-ray photons with energies E 1 > Eth are emitted 
is called the burst duration for the threshold Eth- The 
burst duration as a function of the threshold energy is 




Fig.l. Burst duration versus the threshold energy of 
gamma-ray photons for various evaporation models: (1) 
the model without a chromosphere [4], (2) the DK02 
model with a chromosphere [5], and (3) the H97 model 
with a chromosphere [6]. 
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Fig. 2. Distribution in the parameter Fi for the data 
gathered during nine years; the Gaussian fit is shown 
by the solid line. 

plotted in Fig.l for three evaporation models. Until re- 
cently, the bursts of high-energy 7 rays from the final 
stage of PBH evaporation were sought using several ar- 
rays for detecting the extensive air showers (EASs) from 
cosmic rays [7] and the Whipple Cherenkov telescope 
[5]. Since the threshold energy of the primary 7-ray 
photons is high and the duration of the high-energy 7- 
ray burst predicted by the chromospheric models is too 
short (much less than the dead time of these arrays), 
the results of these experiments can be interpreted only 
within the framework of the evaporation model without 
a chromosphere. 

Within the framework of the chromospheric models, 
the PBH events can be sought directly by the method 
of detecting spikes in the total counting rate of the ar- 
rays operating in the regime of the detection of a single 
cosmic-ray component. The particle arrival direction in 
this regime is not determined, while the effective energy 
of the primary gamma-ray photons depends mainly on 
the altitude of the array above sea level. Such a tech- 
nique was employed earlier to seek cosmic gamma-ray 
bursts with energies exceeding a few GeVs [9j [TUl Qj] ■ 

EXPERIMENT 

The Andyrchy EAS array of the Baksan Neutrino 
Observatory (Institute for Nuclear Research, Russian 
Academy of Sciences) is located at an altitude of 2060m 
above sea level and consists of 37 scintillation detectors. 
The area of each detector is 1 m 2 . To detect a single 
cosmic-ray component, the total counting rate of all of 
the detectors is measured each second. The search for 



7-ray bursts using this technique is carried out against 
the high cosmic-ray background (u> = 11440 s _1 ), which 
requires the highly stable and reliable performance of all 
equipment. The monitoring is realized through simul- 
taneous measurements (with 1-s acquisition rate) of the 
counting rates in the four parts of the array compris- 
ing 10, 9, 9, and 9 detectors, respectively. A detailed 
description of the array and its operating parameters is 
given in [12]. The probabilities P(E 7 , 0) of the detection 
of the secondary particles that are created by the pri- 
mary 7-ray photons with energy E 1 and that incident 
onto the array at the zenith angle 6 in a hypothetical 
infinite-area Andyrchy detector were obtained by sim- 
ulating the electromagnetic cascades in the atmosphere 
and detector [IT]. Since the probability of the detec- 
tion of a 7-ray photon is a fairly smooth function of the 
photon energy, the median energy of the primary 7-ray 
photons detected by the detector depends on their en- 
ergy spectrum. In the case of a source with zero zenith 
angle having a power-law spectrum with an exponent of 
—2.0, the median energy of the primary 7-ray photons 
is equal to 10 GeV. 

Deviations in the array counting rate lasting At < Is 
are sought using the parameter Fi that is equal to the 
deviation (measured in units of the Poisson sigma) of 
the number of counts ki during the ith second of a 15- 
min interval from the average number of counts k during 
this interval: Fi = {ki — k)/Vk. Since variations in the 
cosmic-ray intensity over a time of 15 min are negligible 
in the first approximation and the average counting rate 
is fairly high, one can expect that the parameter Fi has a 
Gaussian distribution with the zero mean value V = 
and unit standard deviation a = 1.0. The parameter 
Di is used to characterize the differences in the count- 
ing rates of the array parts: Di = |X^i=i(-^/ ~ Fi) 2 ■ 
Here, F? is the deviation of the jth part of the array and 
F is the average of the four values for the ith second of 
a 15-min interval. Further processing was done only for 
the 1-s points for which the condition Di < -Dbound (-??:) 
was valid. The values -Dbound(^i) were obtained by the 
Monte Carlo method under the assumption that ki has 
a Poisson distribution. This condition makes it possible 
to reject points for which the differences in the count- 
ing rates of the array parts are unreasonably large, i.e., 
eliminate the instrumental error. The useful events (1-s 
points) can be rejected with a probability of 2 x 10~ 9 
for all of the events and with a probability of 1.3 x 10~ 3 
for events with Fi > 5. About 0.01% of events were re- 
jected from the entire volume of the experimental data 
according to this criterion. 

Figure 2 shows the experimental distribution in the 
parameter Fi plotted using the data gathered during 
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nine years (2005.4 d of the net observation time). The 
experimental data up to ~ 6a are fitted well by a Gaus- 
sian distribution with the mean value V — —(0.0025 ± 
0.0004) and standard deviation a = 1.006. The only 
event with a large (7.9a) deviation was detected April 
17, 2002 at 17:31:29 UT. No 7-ray bursts were detected 
at this time by space-borne instruments. 

EXPECTED SIGNAL FROM PBHS 

Let a PBH be located at distance r from the array 
and observed by this array in the direction specified by 
the zenith angle 9. Then, the mean number of 7-ray 
photons detected by the array is equal to 

^r) = ^N(At,6), (1) 

where N(At, 8) is the total number of the gamma-ray 
photons emitted by the PBH that can be detected by 
the array: 

00 

N(At,6)=[ dE 1 P{E, , 6) ^ ( Ai) . ( 2 ) 


Here, {dN 1 / dE 1 )(At) is the spectrum of gamma-ray 
photons emitted by the PBH during the time interval 
At until the end of the PBH evaporation and S(8) is 
the area of the array. The number of bursts detected 
during the entire observation time T can be written as 

N = p phh TV e g, (3) 

where 

00 

V cS = J dilj drr 2 F(n, n(9, r)) (4) 


is the effective volume of the space surveyed by the ar- 
ray, Ppbh is the number density of the evaporating PBHs, 
and F(n,n) = e~ n n n /n\ is the Poisson probability of 
the detection of n events for the mean value n. 



the number density of evaporating PBHs at the 99% 
confidence level is calculated by the formula 

4.6 

PV±m = TT~^> ( 5 ) 

where the net observation period is T = 5.5yr. The 
substitution of the V c g values for each chromospheric 
evaporation model into Eq. (5) yields a 99%-C.L., with 
the upper limits 1.8 x 10 12 and 1.7 x 10 13 pc~ 3 yr _1 for 
the DK02 and H97 evaporation models, respectively. 

PROBABLE PBH EVENT 

Let us assume that the event of 7.9-cr deviation is 
caused by the 7-ray burst from an evaporating PBH. 
Then, according to the evaporation model [4] without 
a chromosphere, the average number density of evapo- 
rating PBHs is equal to 2.9 x 10 12 pc~ 3 yr _1 , which is 
many orders of magnitude higher than the upper limits 
obtained earlier in several experiments pj [5] ■ 

According to the chromospheric models, the aver- 
age number density of PBHs is equal to p = 7.8 x 10 11 
and 7.2 x 10 12 pc" 3 yr" 1 for the DK02 and H97 mod- 
els, respectively. Note that a single spike in the to- 
tal counting rate was also observed at the EAS-TOP 
array [S]. The deviation amounted to 20.6er for At 
= 2 s. If this spike is a PBH event, then the aver- 
age number density of evaporating PBHs amounts to 
p = 6.2 x 10 12 and 3.7 x 10 13 pc'V -1 for the DK02 
and H97 models, respectively. Despite the significant 
(more than an order of magnitude) difference between 
the mean values, the 99% confidence intervals for each 
model intersect with each other. Thus, if the events 
with large deviations in the total counting rates of the 
EAS-TOP and Andyrchy arrays result from 7-ray bursts 
from evaporating PBHs, then the average number den- 
sity of the evaporating PBHs is 7.5 x 10 10 - 5.2 x 10 12 
and 4.5 x 10 11 - 4.8 x 10 13 pc^vr" 1 for the DK02 and 
H97 models, respectively. 



CONSTRAINTS ON THE PBH NUMBER 
DENSITY FOR THE CHROMOSPHERIC 
MODELS 

The upper limit for the PBH number density (ex- 
cluding the single spike in the counting rate) can be ob- 
tained from the absence of events with deviations > 6a 
for At = Is, i.e., the effective volume of space (4) sur- 
veyed by the array is calculated for n = 6a = 642. If 
the evaporating PBHs are uniformly distributed over 
the local Galactic region, then the upper limit pu m for 
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